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ABSTRACT 

The thermal theory of laminar flame propagation for hyditegen- 
bromine mixtures is described. The method of analysis follows the 
earlier work of von Karman and Milltm and of von Karman and Pen- 
ner. The problem is materially simplified by introducing the steady- 
state approximation of classical chemical kinetics for the concentra- 
tions of H- and Br -atoms. 

i 

The general formulation of the problem is presented in Sec- 
tion 11. 

Approximate solutions of the relevant mathematical problems, 
utilizing various procedures developed by von Karman and his collabo- 
rators, are given in Sections III to V for bromine-rich, stoichiometric, 
and hydrogen-rich mixtures, respectively. 

Numerical values for the physico-chemical parameters, which 
are required for the calculation of absolute values for the burning 
velocities, are considered in Section VI. 

,3k 

Thg calculated burning velocities are compared in Section VII 
with experimental results obtained by R. C. Anderson and his col- 
laborators. Reference to Section VII shows that the variation of the 
observed burning velocities with mixture ratio is predicted correctly 
by the theoretical calculations. Absolute values for the laminar 
burning velocity cannot be estimated with certainty because of the 
possible existence of large errors resulting from extrapolation of 
low-temperature kinetics and heat conductivity data. 
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I. INTRODUCTION 



In recent years a large number of publications have appeared 

on experimental measurements and theoretical studies of laminar flame 

propagation in premixed gases. For a complete bibliography reference 

m ( 7 ) 

should be made to the books of Lewis and von Elbe' ' and of Jost'"' or 

( 3 ) 

to a survey paper by Evans' ' . Experimental measurements on hydrogen- 

M 5\ 

bromine flames have been described by Anderson and his collaborators'*’ . 
As yet incomplete theoretical studies on this flame have been published by 
J. O.Hirschfelder and his collaborators^’^, who are attempting to carry 
out a rigorous numerical solution of the complete flame equations for the 
stoichiometric mixture ratio. 

The method of approach which is followed in the present investiga- 
tion differs in scope and purpose from the type of calculations carried out 
by HirschfeH^r . We follow the work of von Karin an and Vfillan^ and of 
von Karman and Penner' ’ ' in that we shall content ourselves with an 

approximate solution of the flame equations for a thermal theory, on the 
grounds that &he required physico-chemical parameters are not known 
with sufficient accuracy to justify numerical solutions. On the other 
hand, in order to obtain a reasonably severe empirical check of the 
physical model and of the mathematical approximations , complete calcula- 
tions as a function of mixture ratio are carried out. In order to obtain the 
desired results, we shall utilize the steady-state approximation of classical 

reaction kinetics for the chain carriers, as suggested for flames by 

/ / (q) 

von Karman and Penner' , although the steady-state approximation cannot 



be justified for Er atoms at elevated temperatures, it appears that the 



- 2 - 

calculated burning velocities are nevertheless approximately correct* • 
There is a well-known cold boundary-value difficulty whicl 
enters into the solution of the one-dimensional laminar flame equations • 
We shall not discuss this problem except to note that a convenient method 
for avoiding complications is to utilise the concept of an ignition temper- 
ature<“> . 



* This result can be established by an extended analysis camparing 
reaction rates with and without the steady-state approximation. The 
relevant analysis has been carried out by von Kjfrm^n and Millin. 
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II. GENERAL equations for the thermal theory of 

LAMINAR FLAME PROPAGATION IN HYDROGEN-BROMINE 
FLAMES WITH THE STEADY-STATE APPROXIMATION FOR 
H AND Br 



For the theoretical study of laminar flame propagation in hydrogen 
bromine flames, we shall utilize the accepted reaction scheme for isother 
mal production of HBr in stationary systems at temperatures of the order 
of 500°K. The required equations ard:^*^ 



X’ + Br 2 


2Br 


+ X’ , 


(1) 


Br + H 2 JX 


HBr 


+ H , 


(2) 


H + Br 2 


HBr 


+ Br , 


(3) 


H + HBr^ 


H 2 


+ Br , 


(4) 


2Br + X • 


Br, 


+ X 1 


(5) 



Here X' denotas the "third body", i.e. any of the chemical species 

rk 

H, Br , or The specific reaction rate constants to 

kg for the five reaction steps are indicated above each of the arrows. 

It is well-known that the steady-state approximation for H- and 

( 12 ) 

Br-atoms leads to the relation' ' 



d(HBr) 

3t 



2k 2<E 



1/2 

) (H 2 ) (£r 2 ) 



(HBr 
k 3 '^ r 2 



1/2 



( 6 ) 
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where the chemical symbols in parentheses denote concentrations for 
the cherr ical species and t represents the time . For thermal flame 
theories it is convenient to introduce mole fractions, X; we identify 
H3r , Br^i and by the subscripts i, 2, and 3, respectively. Then 

(h 2 * ’ifr x 3 *#r (x 3,o < 7 > 

and 

(Br 2 ) X 2 = A (X 2,o “ “F> (8) 



v/here p denotes pressure, R is the molar gas constant, T is the temp- 
erature, and the subscript o identifies initial conditions. 

For flames without inert diluents it is convenient to identify the 
initial composition by the parameter a. Let X, = a, X, = 1-a. Then 
Eq. ( 6 ) becomes 



P d < X l 

Ft “3T: 



1/2 3/2 x x 3/2 

) Zk 2 { iF ) (a • -r ){l -~~r ] 



(9) 



kj) x i 



If Wj , an{ * ^ 3 ^ eno ^ e the molecular weights of HBr , Br^ . 

and . respectively, then 

d(x i> » 

v 'iTT ivr 



represents the rate of change of the mass of bromine per unit volume 
with time. But this quantity ir> also equal to 

w i d h 

m — — 

dx 

if in iG the mass flow rate per unit area of the mixture (i«e 0 the eigen- 
value for the problem of laminar flame propagation), x equals distance, 



and 
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vV =Z W.X. 

1 i l 



( 10 ) 



is the average molecular weight of the gas mixture. Hence Sq. (9) 



beco neG 



1/2 



3/2 



X, 



dX, 
m 1 



(a- -=r-)(l -a- -^4 

d. a) .(^.^l) Xl 



X 3/2 



( 11 ) 



In the absence of diffusion the general energy equation for one- 

jo) 

dimensional flow can be written in the form' 7 



^ * (9 - D 

me 

P 



(Y. - Y. ,)h° 

+ -i-lj hf -L 



( 12 ) 



c T. 
P * 



where A equals the thermal conductivity^ the reduced temperature is 
0 = T/T^, with the subscript^ identifying equilibrium conditions at the 
adiabatic flame temperature; Y^ is the weight fraction of the j'th chemical 
species; h . equals the standard specific enthalpy of formation at temper- 

J 

.at. 

ature T q of the j'th species; the average specific lieat"c^ is given by the 

• (9P~ 

expression ' ' 



s f <Y j.< ‘ Y j.o ,h ° • 



03) 



For the present approximate calculations we shall ignore the 
production of bromine and hydrogen atoms at the equilibrium conditions.* 
In this case it is found from F.qs. (12) and (13) that 



* This assumption is in line with the steady-state postulate but is not 
correct at elevated temperatures where appreciable concentrations of 
bromine atoms may form in bromine-rich mixtures. The required 
corrections have been worked out and have been found to be relatively 
small . 



since h^ s 0* But 



< x 2- x 2.f> = * x r x M> 



and 



< X 2.f • X 2.o» = Z< X l,f- X l.o> 



Hence 



A 

m ' P 



d6 

355 



( W.hj + 1 W,h°) 

(6^1)-(1.9 o ){X r X 1>£ ) i I 7 2 w7^) ' 



We shall be interested only in flames for which X, = 0. Hence 

1, o 

the energy equation reduces to the following simple expression: • 

X 

•4- ® “<»-»■ • < u > 

me? l,f 

Following the conventional procedure for solving problems in 
laminar flame propagation, we now divide the continuity equation 
given in Eq. (11) by the energy equation given in Eq. (14) in order 
to eliminate the distance x as the independent variable. In this 
manner it is found that 



- 7 - 



m 2 ^ dX 



? *35“ = 2k 2 ( F7 } 



x 1/2 3/2 



3/2 



AW UB * "5 

We now introduce the assumptions'' 7 that 



jl-a-4- ^|)X 1 ][(9-1)+(1-0 o )(1-X 1 /X 1>£ )| 



1/2 

A - A f 9 


(16) 


1/2 -9 ? /9 

k 2 =B 2 9 e , 


(17) 


kj B £ -Qj/0 

77” =s— e 
*'5 B 5 


(18) 



Equation (15) now becomes 



dX £ 

*35" 



A f W B. 1/2 
■ST’ (g-> 9 



3/2 -1/2 -9’ Z /Q 



m c 



(19) 



Ua-X^l-a -^) 372 

X 



|(1 .al,( ‘ . ^xj [(O-lW-^U-xi )] 



wnere 



«' = 9„ + ± 9, 



( 20 ) 



'2~' , 2 r 7 w l * 

It is now convenient to introduce the parameter -/j. through the 



relation 



B, 



1/2 



A f W "l 

^ "7W B2 ^5 ) t frf ) 

me 5 i 

P 



3/2 



• 05) 



( 21 ) 



whence Eq. (19) becomes 



3/2 



dX 



1 Ji 0 



-i/2 -0^/0 



■ar 






(e-i)Ki-e )(i-^-) 

° x i,f 



(l-a)-4 - M k 3 )X l 



( 22 ) 



The boundary conditions are: 

9 = (h when X ^ = 0 



(23) 



and 



0 = 1 



when 






= 2a 



= 2(1 



for a < l/2 
for a = l/2 
a) for a /> 1/2 



>• (24) 



The solution of the boundary value problems is discussed in the 
following Sections III, IV, and V for a < 1/2, a —1/2, and a y l/Z 
respectively. 
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BROMINE -RICH FLAMES 

For the bromine-rich flames, solution of the boundary-value 
problem can be effected conveniently by use of the von Karman- 
Millan technique. ^ The boundary conditions are 



and 



Q = when Xj = 0 



0 = 1 when Xj = 2a 



(24) 



From Eq. (22) 



1 



■ dx. , T - 9 2 /9 - A 3/ 2 

W ~ 9 .. !■ Jt - t 2 *-**)* 1 -*- 2 ) 



(9-I).(l-S o )(l-|l) (l-a-OJSlXj) 



(2S) 



where the appropriate numerical value of k^/k^ has been introduced 
Use Section VIC for details). 

Rearranging Eq. (25), we get 



„„4 -»2« d« O-a.-0. ssxyp 

6 V-'a^X,) 5 ' 2 



X, 



(e-DHi-yd-^) 

2a-X . 



(26) 



from which it follows that 



-O' 
w 2 



,, * „ , d0 x (1-1. 762a) LO-4-1 

Tie lim(_)=l ■ - L T7f ~ — ~ 

9 — <► 1 ^1 (l-2a) J/<: 



lim (4£-) - (1^) 



^ 1 ~ 



-i 



(27) 



or 



lim , d0 4 



Ae’’ 2 

(l-2a)' 3/i 



<1 ~°o ) (l-1.762a) 



"2a" 



(1 "2a) 



Tff. 



(28) 



Defining 



i t 1 -1 . 762a 

6 (a) YfT 

(l-2ar^ 



( 29 ) 



-10- 



and 



1-3 



b ' < a - V ° -25— 1 



(30) 



then Eq. (28) becomes 



lira , d8 , 6(a) b (a, 3 o ) 

e->l ( 3!T )= ^5 

1 Jle c t<5( a) 



(31) 



Introduction of Eq. (31) on the left hand side of Eq. (26) for 



8~1, 2a, leads to the following relation: 



J te 



-3’ 

y 2 



X, 



b^a', 6 q ) 

Tie 2 + 6(a) 



(e-D+u-s^d-^i) 



2a-X, 



or, using Eq. (30), 

1-3 = b'(a,8 o )(2a-X 1 ) 



1 - 



Tte 



>2 



— ^sn; 

Tie 2 rfi(a) 



or 



1-9 '= b'(a,8 o )(2a-X 1 ) 



6(a) 



Iff 1 " 

Tie 2 + 6(a) 



(32) 



Introduction of Eq. (32) on the right-hand side of Eq. (26) 



yields 



A 3 



\ ' 9 2 /8 d8 






(1 -a-0. 331Xj) 

. _r--57r 

(l ~ a ~g Xj) 



b’(a,9 o ) 



6(a) 

““^2 
Tie 6(a 



;(33) 



integration of this last relation between © ft , and 9 = 1, = 2a, 

2a 

_^9 0 )e ” 2 | 



leads to the result 

•1 1 



$ 



B -Z e - a 2 /B d8 a b'[^» 0 )e 

c Ae 2 +6(a) 

y t 



X 



(1 -a-0. 35AX.) 

/• a T x An dX i- < 34 » 

(l-a-^ Xj) 



t 



- 11 - 



For 9/^C^ 0, we get from Eq. (26) 



■A. 0 



4 - e 2 /9 



d9 



<e-i)+(i-e o ) 



2a 



(l-a-0. 331 Xj ) 

ft“37f 



-* Xj) 



(35) 



the integral of Eq. (35) from 0. , 0 to 0^, X^ is 
>0. 



* -7 ‘ 8 2 /8 

-2 2 de = ■y——r 

2a vt 



Qi (9-DH i-« 0 > 



(l-a-^Xj) 3 ' 2 1 



(36) 



But 



( 8 ) 



3 A ^ Cl /A 

9 * e 



9 i 



i -I Vi (1 -ST de 






r e t i 

J 8 '^' 



8 2 /8 



d@ 



whence 
•0 



J 






- 8 2 /8 , (# t - IW »-» 0 > 

c de ~ — 

2a A. 



"X. 



(l-a-0. 381X.) 

— j -g-rydX.. (37) 

o < j - a 4 x i> 



d0 

Sincere— is to be continuous at 0.»X. it follows from Eqs. 

Q v w 

(35) and (33) that 



<0 t -lHl-9 o ) 

zs 



j (l-a-0. 381%^ 

1 

e 



nr37r 



A.Q 7e 2 fc ^’ a ‘I X l ) 



l-a-0. 381X 



1 



n ' 1 y. .3/2 -1 -6«/6V 

(I-a-^Xj) A9 2 



b'(a,0 Q ) 



1- 



iial 



Tie ^i^(a) 



whence 



<0 t -l)+(l-9 o ) b'(a.9 o ) e 



-0' 

y 2 



"®2 

Tie + S(a) 



2a vl 



( 38 ) 
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In view of Eq. (38), Eqs. (34) and (37) can be added to yield 



the result 

1 



_q» _ 2a 

-i -91/0 b'(a, 9 ) e 2 \ ( 1 -a-0. 381X.) 

i ^ _ O 1 A i v 

j rz ’JTz dx i * 



9 “ e de = 



-IF 






2 + S(a) 0 <1 * a ”I X l)‘ 



In order to evaluate the integral 

• 1 



S 

9. 

i 



-j -«2 /» 

# e d6 



we set 



1 12 1 4 . 2 . 

= y, 3 - = y , —J- = y and d9 = - dy 

W " © y J 



Substitution of these relations yields 



j 



7 



9 e 



d9 = -2 



But 



/ 



1*1 



-o;y 



1 

•r 



y c 



4 dy = 2 

y 



-9'V 

e ^ dy. 



*T 



c/ 



-Q'y 2 

e dy = 



*7 



■EOj ye 



-9'y 2 






-2y0' 



or 



-9'y 1 



2y9* 



e 






j i 



, v 

y 

1 



^1, -e.y 2 

~~7T e dy» 



( 39 ) 



Comparison of the two preceding relations leads to the result 




(8) 

For sufficiently large values' ' of 6!, /Qj we obtain* 



1 






1 * -9‘/9 / -9L I I t 

0 2 e 2 dO = ( 2e 2 -2 *57 e 2 J -2 [2 \e" t dt 



0 0 



r -t ; 

T 



= ^2g" 62 - 2 >57 -2>?5£ I - E 2 CITST )1 



where 



<7>i - 



y* 



2^2)= T E 2l' w 2 



f -t 2 
= \ e dt. 



But 






-1 


-9' 

2 




^ e 


ZW Z 



, 1 , 1.3 1.3.5 

1 " Tip' + ~z 5 + 

Z ^2 (ZQg) <2e^r 



*In the practical applications considered here, Q^/Oj * s su ^^ c ^ en ^Y 
large to justify the approximations made here. 



whence, finally, 



-l -ej ,/e -e* -o' 

0 c e d0^2e -2e 



0. 

l 



1 , 1-3 1 • 3* 5 , 

JqT + ? g-+ . . . 

C *Z (20!,r <29^r 



or 
'1 



-i -01/9 ’®2 

S e dO - -4n — 



, 1-3 L 1-3*5 

W"' 



(40) 



- 9 2 



vie 



-91 



Equation (39) can now be solved for vie with the result 

»-• ■>.''* r 2a U-a-0- 381X1)^1 



(1-a-^Xj) 



IIP r 

W 2 r 

2a — 
w 2 



, 1-3 L 1-3-5 

‘■i 25 ? W" 



1-1. 762 a 
(l-2a) S/2 



or 



71 e 



) (l-a-0. 381 X, ) 

0 n"a-i/^^p^ dx i 



2a 



1-3 

'(2^ 



1 ^ r ri k t 

Z ] (20‘) 



1* 3* 5 



(1-1. 762a) 
(l-2a)*^ 



(41) 



where <5 (a) and b'(a,9 Q ) have been replaced according to Eqo. (29) 
and (30) respectively. It is now necessary to evaluate the integral 



.2a 

r. d-a- 



(l-a-0. 3813^) 



1“'T37I 

1 x i' 



dXj = 



s* ucl 



-2a 



ii-a* 

T*l> 



T7I 



dX x - 0. 381 



Xj dXj 






E\ r idently 



-lb- 



,2a 



(l-a)dX. 

" 1 . .3/2 = 4 ^ l_a) l 

(1-a- jXj) 



1 



nrzr rrz 



and 



0.381 



2a 

5 



XjdXj 



nrvs 7 i 



= 0. 381 



(1-a-^Xj) 



8a 

yra- 



ie(y 1 -2a -yr< 



Therefore 
2a 



^ (l-a-0. 



381X.) . 

"TTY dX, = - • 

(i-a-^Xjr 7 ^ yiT2F 



4(1 -a)~3. 05a-6. 10H2 



. 20aj 



yr^ 



M 



4(1 -a)-6. 10(1 -a) 



= 5 .* l 5a * 2 ‘ i£ + 2. 10(l-a) £ 
X2 -2a yi-a 



The final result for bromine -rich flames becomes 



A. e 



-O' 

¥ z 



(i-e o )oy(a) 



2a 



, 1-3 . 1-3*5 

l ~(7on ■ j 

KC *Z' (28 ' 2 r 



(1-1. 762a) 
(l-2a) 3/ * 



(43) 



from which the eigenvalue A can be obtained for bromine-rich mix- 
tures . The mass burning rate m is then determined through 7 "q. (21). 



V 
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I\'. THE STOICHIOMETRIC MIXTURE 

£ or the case a = •j and X^ ^ = 1 the basic differential equation, 
Eq. (22), becomes 



* -0* /9 

dX l -A/yi9 I c 2 

*ar = 



5/2 

(1-Xj) 

irromqr 



or 







(1-Xj) 5/2 

(1-0. VbiXj) 



d9 

dX x 



= (0-iHi-e^i-Xj) 



(44) 



In order to integrate Eq. (44) we apply an ingenious, as yet un- 

/ / 

published, method invented by Dr. Theodore von Karman. Ec replace 
the singular point (0 = 1, - i) by (9 = 1, X^ = Xj< 1); then Eq. (44) 

becomes, for (9 - 1) small compared to (l-9 Q )(l-Xj), 



-7 - 9 > /Q l o 

A 9 2 e 2 d9 1 '®o 

T=wn snq" 



(45) 



Integrating Eq. (45) from (9* , X ? ) to (l,xp leads to the result 

f 

i t 

ft I rs - w j (i-Xj)- 



-9^/9 

d9 = (I-9J 

yir ° 



,X{ 



(1-0.761 X.) 
T7I dx i 



9* 



V 



(46) 



For 1, the preceding expression becomes 

i 

l 

(1-0.761 X.) 

t. 



x; 



_c 

~e 2 (1-9 *)o>( 1-9 o ) 



X f 



(1-Xj) 



TTT 



dXj . 



(46a) 



If Eq. (45) holds, a© a first approximation, from (9^,0) to 



(i.XJ). then 
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f -®J/® 

— 1- de = (l -o ) 

y?j w 

e. 0 



X 1 



( 1 -0. 76lXj) 
d-Xj) 



7^71 dx i 



(47) 



or, utilizing the result obtained previously for the integral appearing 
on the left side of Eq. (47), 



-9! 



A. e 



y 1*3 1*3*5 

f* 5 ? 






x i 

\ (i-a76ix 1 )dx 1 

J (i _ Xl ) i/2 



(48) 



From the original differential equation it is apparent that the 
limiting slope at the hot boundary io 



H: ° 



to which we join the curve of Eq. (46a) at (9*,X^). Therefore, at 
(Q*,Xf), 



(l-0*)»(l-0 o )(l-X*) . 



(19) 



From 2-qs. (46e) and (49) 



- 9 ! 



— e 2 (l-X.») = 

n 1 



y» 
X 1 



X l* 



(1-0. 76lXj) 

d-Xj) 



TTT dX i 



r .^ de 



With = (l-9 o ), Eq. (45) becomes at (9*. X^): 



(5f)) 



vt e 



-9! 



3/2 



ti- 6\ vsref) {l ~ x V 



= 1. 



( 51 ) 



7/e next integrate Eq. (50) noting that Xj* and XJ arc very 



close together so that we may approximate (1 - 0. 761X^) in the inter 

val Xj£ to XJ by (1 - 0.761X^). Hence Eq. (50) becomes 

X* r 

^i-x^ -y/nx[ 



•0J 



~ e 2 U-X«|) = (1-0.761X1) 



f l dX t 

^ (i-x,) 

XI 1 



■ S J 2 = 2(1-0.761X1) 



/(1-X*)(1-X1) 



or 



A 



"°2 (l-Xf ) 3/2 y/TTxf -y/TSJ- 



* 2 1 



KToq* 



(52) 



From Eqs. (51) and (52) we now see that 



2 F*T 2 _, 

2 /l^q" ' 2 - 1 



or 



/ 



Tocf 

t^t 



3 

I 



whence 



x ! =! x i-| 



(53) 



The next step is to obtain an explicit expression forXj from 
Eqs. (48) and (50). Equation (50) may be written in the form 



A ± 
n 



•O' 

y 2 



X 



1 C (1-0.761 X 

(iTxT^- 

XI 



Xi)dXi 



(54) 



or, using Eq. (48), v , 

(1.0 o ) > X MI-0.761X 1 )dX l 



ILiM 



372" 



f , 1-3 , 1-3-5 

I 1 -/ tkt : + . : ' 2 



P^i-Xf ) 



(l-0.761X 1 )dX 1 

rs~r^ 

(1-X^^ 



Replacing X^ in terms of XJ according to Eq. (53) in the preceding 
expression leads to the following explicit relation for XJ: 



, 1*3 1-3-5 

1- — 



<2°2 )‘ 

|i*- x i> e 2 



(1-0. 761Xj) 



9 x* 5 

? X 1 ? 



U-Xj) 



37T dx rd*V 



(1-0.761:^) 
(1-Xj) 



T7F dX i 



(55) 



Integration of Eq. (55) leads to the following first approximation 



for X{: 



X i = l 



f. 1*3 ,1*3*5 1 

liSP"": 1 



V (29^) 2 
T- 6 2 <1 ” 9 o ) 



(56) 



Hence, from Eq. (53), 
X^ 1 



and, from Eq. (49), 



^0*cri - 




(57) 



(58) 



Comparison of Eqs. (56) and (57) shows that Xj and X^ are indeed 
so close together as to justify the approximation made in the deriva- 
tion of Eq. (52). 

' 0 2 

In order to obtain a first approximation to Vt e we proceed 
from Eq. (54) as follows: 



^ -9£ 2(1-0. 76lXj*) 


f i . » _ 


ft C 1_X t 


t/roq- yrxf 


1-0.761X^ 


l-0.761(^X' - |) 


- "■ 37 1 ~ ' 

(1-X-J) 


r 9 ~5 t 3/2 

[1-(| Xj - | ) j 



after replacing Xj according to Eq. (56) it is found that 



A e - e 2 

n 



0.239+0.761 



1- ill ^ kJi ■> 

‘(29!,) (io^j 2 

k^oT 



f. 1-3 .1-3-5 1 

rcg? W" J 



T7z * 



(59) 



30 2 (l-0 o ) 

Therefore, for the i*ange of values of 0i, actually encountered, 

-9' 

Ae )T5 yS^i-Oo) [ 0.761 + 0.717 0^(1-© Q > J (59a) 

from which the eigenvalue A for the stoichiometric mixture ratio 



can be calculated. 



V. HYDROGEN-RICH MIXTURES 

The mathematical technique which we have employed for 
stoichiometric mixtures can be used also to solve the boundary- 
value problem for hydrogen-rich mixtures. However, since the 
calculation is somewhat laborious, we follow another suggestion 
by Dr. von Kifrm£n, according to which v/e test a more highly ap- 
proximate procedure for mixtures containing 43. 3°/o Brg, for 
which the solutions become particularly simple. The success of 
the simplified procedure for mixtures containing 43. 3°/o B^ gives 
us confidence in the results obtained for hydrogen-rich mixtures, 
in general, for which v/e utilize a method of calculation in which 
the limiting slope of the C vs. Xj curve does not have the correct 
value at the hot boundary. 

A. The Mixture Containing 43. 3°/o of Bromine 



simplification is noted in the basic differential equation. For 
these mixtures X. - = 2(1 -a) and it is seen that 

~ I 



Equation (22) becomes, for (9-1) small compared to (l-9 o )(l-Xj), 



For mixtures containing 43. 3°/o of Br^ a considerable 



•, 1 -a 

2 ** 03ffT • 



1 -a = 0.433; 




X, 3/2 




- 1/2 - 9^/9 



( 60 ) 



e 



(1) Direct Integration of Eq. (60) 

■Without regard to the lirpiting slope at the hot boundary, 



we can integrate Eq. (60) directly with the result 



J l 



1 

i 






2(l-a) 



-e*/Q 



dO = 



0.381(l-0 o ) 

Z7Z5T 






dX, 



27171 

(1-a-^i) 



( 61 ) 



since 1-a = 0.433. Hence 



, 1-3 . 1* 3* 5 

lm &2 + <15y*’- 
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Tie 
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i j;3 1*3*5 _ 


0. 381(l-0 o )(-4) 
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e* 
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o 
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- 2.318(l-0 o ) 




and 


-0» 2.318(1 

-A e a — 


• 

-1N3 

O 

CD 

1 





2(l-a) 



( 62 ) 



The fact that*Eq. (62) yields a fair approximation to the solution 
of the boundary- value problem will now be established by a more 

it 

careful integration in which d0/dXj has its correct limiting value 
at the hot boundary. On the basis of this result, we shall inte- 
grate the flame equation for all hydrogen-rich mixtures directly 
in Section VB. 

(2) Integration of Eg, (60) With the Correct .Limiting 
Slope at the Hot Boundary 

Integration of Eq. (60) from (0*, XI ) to (1, XJ) leads to 



the relation: 
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- 0 2 /e r 1 dx 

l dC = j " 1 " a ”'^i.r/2 • 



0* 






(63) 



For 6* ^ 1 



-9' 

Tie 2 (1-0*) = 0.8S0(l-6 o ) 



dX, 



5TT7I 

Xf (1 ' a ” T ] 

If Eq. (63) holds from (0^,0) to (1, XJ), then 

l X! 



A 



e 



-O’ 
*2 



T TT 



d0 



C‘ 

= 0.880 (1-0 ) J - 



dX, 



0 4 0 



Xj 1/2 



0 (1-a- ) 



(64) 



or 



X* 



A 



> 0 * 

*2 



Q* 

w 2 



, 1*3 , 1*3*5 

<**P w;?"" 



? 0. 880(1 -0 Q ) 



J 



dX 



1 



“ XTT72 

(l-a--j~) 



(65) 



We now join Kq. (64) to 

-* 1-8 

d0 1 *o 

a3?j 



at (®$. X>|) , 



since dO/dX^ = (!<-0 o )/2(l-a) is the correct limiting slope obtained 
from Eq. (22). At 0 = 9*, Xj = X>!(: 



1 - 0 * 

zcr-aj - ^ 



1_6 o 



(66) 



and, from Eqs. (64) and (65), 
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"fl-ai" = 0.880 ^ 



dX, 



X 1 



r ^IF 2 



With 



d0^ 

cnq- Krsr 

Eq. (60) becomes at (6*. X*|) for 0* 

X* 1/2 

X 1 



i: 



Jle 



(1-a-y) 



= 0.380 



(67) 



(63) 



From Eq. (67) 

JLe 



“®2 [2(l-a)-X*] 

— ztt-T)— 



= -3. 520 < 






*, 1/2 



(69) 



But, from Eqs. (68) and (69): 

i 1/2 

fZ x 0. 880 [2(l*a)-Xjf ] = 3. 520 



t x r i i/2 


r X'l 


1/2 


[<*-»- J. ] 
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■ 1 (1-a)- -j- 


> 

4 
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Xf 
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1/2 



.,, 1/2 



( 1 -a)- 
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*, 1/2 



whence 



< l — $> 4 



(70) 



From Eq. (67): 
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‘°2 



X i 



Act * _ 2(1 -a)0. 880 
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X7TT/2 

X{ l (1 ' a, ‘Tl 



and, using Eq. (65), 



0.880 (1-S 0 ) 
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T.T"~F‘3*T 

r* 5 ? <1^7 





Thus ^ 



(l-a> 



, 1-3 , 1* 3* 5 

l *<»? ii^F 



O-V ®2 



R 71 
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x f ^1 
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1/2 



x i 



1/3 1 



(1 -a) - (1-a- -p~) 



Using Eq. (70) to eliminate Xj, the preceding expression becomes 



. f. 1*3 ,1*3*5 ] 

| (1-a) -<3Bp (IgT?' ••• j 
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which yields 

X* 1/2 1/2 

(l- a > ) s (1-a) 
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1 



1.3 , 1*3-5 

(^6*) TtqTT^ 



’2' (29^) 



1 / 2 '| 



Combining Eqs. (68) and (71), and noting that (1-a) = 0.433, we get 



I 



(71) 



_/le 



■0* 



1.16 



1 - 
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. 1*3 ,1.3*5 



L/2 



(72) 
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The fact that Eqs. (62) and (72) yield simila .* * results can 

be seen by expanding the denominator of Eq. (72) for large values 

of (1-0 ) 01. In this manner it is found that 
’ o Z 

Tie c Ci 2. 32(l-0 o )0^ (72a) 

which is seen to be practically identical with Eq. (62) for large 
values of 01,. On the basis of this result we are encouraged to 
integrate the basic differential equation for hydrogen-rich mix- 
tures directly between ©j, 0 and 1,2(1 -a) after neglecting the 
term in (0-1). Details concerning this work are described in 
the following Section VB. 

V.B. Hydrogen-Rich Flames with Arbitrary Concentrations of 
Bromine 

The jre suit for the general case of hydrogen-rich flames 

is obtained moat conveniently by using direct integration of the 

* 

differential equation with the term in (0-1) neglected. For 
hydrogen-rich flames a>^- ans * ^l f = 2(1 -a). Equation (25) 
becomes 
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|l/2 



_ [ 1 - 1 . 7 ( 
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cos 
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* 2 YT-1E 



and Eq. (73) becomes 
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yi x 0.381 -Q rgf 



l-l. 762a 2 co8 -l JSr + 2 /Z=Za 



TOST 



y3H72 



(74) 



Since 0 Q is not a sensitive function of a we shall treat it as 



a constant. Then 



-®2 /e 

dO = 0. 539(1 -OJ g(a) 

o fS ° 

where g(a) is defined as 





1-1. 762a 



2 

/la^T 



cos 



■Wl 



+ 2 



(74a) 



(75) 



Finally, integration of Eq. (74a) yields the result 



./L e 



-S’ 

W 2 



°. 539 (!-9 f | 9, g(a) 
_ i»3 1- 3* i> 

r (i^ 



( 76 ) 



from which the eigenvalue -4_ for hydrogen-rich mixtures can be 
calculated. For (1-a) = 0.433, 1.762a = 1, and Eqs. (75) and (76) 
reduce, of course, to Eq. (62). 

'°2 

We have now obtained analytic expressions for Ae for 
all possible hydrogen-bromine mixtures. Before summarizing the 
basic expressions for the laminar burning velocity, we proceed to 
consider the relevant physico-chemical parameters which must 
be used in order to obtain numerical results. 
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VI. PHYSICO-CHEMICAL. PARAMETERS 

A, Equilibrium Constant for the Reaction Br> ^ 2Br 

The value o£ the equilibrium constant K p (i. e. , of the equilib- 
rium constant expressed in terms of partial pressure ratios) for the 

bromine molecule -bromine atom system have been taken from the 

(13) 

Bureau of Standards compilation' ' which covers the temperature 
range between 298°K and 1500°K. Equilibrium constants at 1600°K 
and 1700°K were obtained by extrapolation of the Bureau of Standards 
data. Equilibrium constants as a function of temperature are sum- 
marized in Table I. 

The data given in Table 2 have been fitted to an expression 
of the form 

K p = BTe’ A/RT 

by plotting log(K / T) vs. (1/T). The plot is shown in Fig. 1. Ap- 
P 

propriate numerical values derived from Fig. 1 are 

A = 45, 210 cal. /mole and B = 230. 5 atmos. /°K . 

.a 

The equilibrium constant K , expressed in terms of concen- 
tration ratios, for the reaction 

Br 2 =!=^2Br, 

is related to K p through the expression K c = K p /RT. 

Hence it follows that 

^ _ 230. 5T .-45,210/RT 
K c - SZJ5TT e 

= 2.803 e” 45,2l0/RT moles/cm 3 

3 

if concentrations are expressed in moles /cm . But 
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where kj and kg arc the specific reaction rate constants for the pro- 
cesses described in Eqe. (1) and (5), respectively. 

In terms of the parameters needed for the calculation of lami- 
nar burning velocities we find, therefore, that 

3 

Bj/Bg = 2.808 moles/cm (77) 

and 

Q x = 45, 210/RT f . (78) 

B. The Specific Reaction-Rate Constant for the Reaction Between 
Br and I 4 2 • 

The specific reaction-rate constant for the process 

^2 

Br + H 2 — =-* , HBr + H 

is based on careful experimental studies carried out by Eodenstein 
(14) 

and Lind. ' ' The experimental results have been re-evaluated by 

#12) 

Pease' ' who finds that 

' k 2 = 4. 56 x 10 l ^Te" 18, 7U0/ET cm 3 /(mole-sec), 

(79) 

-it 

whence 

B 2 = 4. 56 x 10*^ )'T f " cm 3 /(mole-sec). 

We shall utilize Eq. (79) for the specific reaction-i ate constant 
in preference to slightly different forms proposed by other investi- 
gators.^’ 7 ^ However, it is important to note that alternate correla- 
tions of the available experimental data lead to greatly different pre- 
dicted values for k 2 at elevated temperatures. Thus Hirschfelder 
and his collaborators^ recommend the relation 

k ? = 3.45 x 10 10 T e" l6 ’ 640/RT cm 3 /(mole-oec). (79a) 
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A comparison of the calculated values of k^ determined from Eqs. 
(79) and (79a) is shown in Table II for the temperature range which 
is of interest in connection with calculations of laminar burning ve- 
locities for mixtures. Reference to Table II shows, for 

example, that the values of k 2 calculated from Eq. (79) at 1500°K 
could be too large by roughly a factor of 1.7. Hence the laminar 
burning velocities, which vary as the square root of k^, could be 
too large by perhaps a factor of 1.3. 

C. Ratio of the Reaction Rates k^/k^ for the Removal of Hy- 
drogen Atoms 

The ratio oi the reaction rates k^/k^ for the removal of H- 

atoms according to the processes 

H >Br 2 HBr + Br 

and k^ 

H+HEr— UH Z + Br 

( 1 ^-) 

is again based on the experimental work of Bodenstcin and Lind' 
and Bodenstein and Jung. ^ ^ The value 

kjj/kj 3 1/8.4 (80) 

does not vary with temperature over the relatively narrow temper- 
ature range for which experimental measurements are made. We 
assume that this result can be extrapolated from about 600°K to 
temperatures as high as 1700°K. Since the activation energies for 
the reactions (3) and (4) are probably very small, this extrapolation 
should not introduce a large error even if the activation energies 
for reactions (3) and (4) are not identical. 
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D. The Adiabatic Flame Temperatures and Equilibrium Compool- 

tions for fixtures 

Adiabatic flame temperatures and equilibrium compositions 
have been determined by tlie use of standard procedures. The 

results for T f as a function of composition of the gas mixture are 
plotted in Fig. 2. ; the equilibrium mole fractions of H T r, H,, Brj, 
and Br as a function of composition are shown in Fig. 3. Relevant 
numerical values are summarized in Table IQ. 

E. Thermal Conductivities In Mixtures 

Values of A/’ as a function of composition for mix- 

tures have been calculated by Cooley and Anderson^ on the basis 
of the elementary kinetic theory of gases, which does not include 
proper allowance for the influence on thermal conductivity of energy 
transfer between internal and translational degrees of freedom. The 
results of Cooley and Anderson, which represent lower limits for 
the thermal conductivities, are summarized in Fig. 4, where we 
have plotted as a function of composition. Reference to Fig. 4 
shows that increases slightly with T£ and more rapidly with the 
concentration of in the gas mixture. 

Actually the problem of estimating appropriate thermal con- 
duction coefficients is greatly complicated for flames by the inevi- 
table composition changes occurring during reaction. Campbell and 
( 6 ) 

Hirschfelder' ' propose a relation of the form 

A = A 1 +A| "T, (si) 

where both and rl" are parameters which must be determined 
separately for each gas mixture. Our formulation of the problem 
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of flame propagation doe3 not permit the use of this two-parameter 
expression for A without a complete recalculation of the results, 
because of the change in the assumed temperature dependence of A. 
In viey/ of the crudeness of the approximations for A , we do not 
feel that the additional labor is warranted. For this reason we have 
used the following relation^ ^ for the thermal conductivity of a 
mixture of species i and j: 



'^f mixture* x i* ~ 1 



AiVAjX. f 



i^j 






(82) 



Equation {82) yields reasonable results for binary mixtures in which 
the pure components have greatly different thermal conductivities. 
Thus Eq. (82) i3 particularly well suited for making estimates of A 
in hydrogen-rich flames. In order to calculate the thermal conduc- 
tivities of the pure components, we utilise the following relations 

( 6 ) 

proposed by Campbell and Hirschfelder, ' ' and neglect again the con- 
tributions of A Br and H at T £ to the thermal conductivity: 

(S3) 

^|Br 2 (T £ ) s 1.20 x 10" 5 , (84) 

6. 08 x IQ" 6 ff 7 

/VHBr(T £ ) ^ (85) 

, A 1.275 x 10 J 
1 + — 

where ^ . is expressed in cal/(cm sec °K). 

Numerical values for A £ as a function of mixture ratio are given in 
Table IV and are plotted in Fig. 4 as a function of composition. 

It should be noted that our method for calculating <A £ 



AH 2 (T £ ) = 5. 13 x 10 



-4 
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cmphasizes the numerical values at temperatures close to T^. This 
procedure receives partial support through the well-known result that 
most of the chemical reactions in a flame occur at temperatures 
close to Tj. 

F. The Average Specific Heat 

A more exact expression for the average specific heat than was 
used in the derivation of Eq. (14) can be obt lined from Eq. (13), viz. , 



-1 



Wl. 



w- 



c = — ’X, ,hr < — (X, ,-X, )h, 

p A £ (i-e o ) L W 1 * f 1 w 2,f 2,0 2 



W 4 o w 5 ol 

w 4, f 4 5, f 5 J 



(85) 



since X, = X. = X c =0 and h, = 0 if the subscripts 4 and 5 iden- 
tify, respectively, the species H and Br. Although the contributions 
of H and Br are small, they have been included in Eq. (85). Average 
specific heats as a function of composition have been calculated both 
with and without allowance for dissociation into hydrogen and bromine 
atoms. The results are listed in Table V. Reference to the data given 

J* 

in Table V shows that the influence of the atomic species on the aver- 

% 

age specific heat is negligibly small except for near-stoichiometric 
mixtures. 

If dissociation into atoms is neglected, then simple explicit re- 
lations can be obtained for the average specific heats. After introduc- 
tion of appropriate numerical values* the following relations are ob- 
tained: 



for Xj £ = 2a, ^2 { ~ ^Z o = O -2a)-(l -a) = -a, i. e. , for a <: ^ 

,4 



2.466 x 10 a 

p' (TO':a-r5-?.Sa)T f (r-& o ) ; 



c = 



( 86 ) 



♦ Appropriate numerical values are 

Wj h ° = -8660 cal/mole and h° = 7340 cai/mole 



( 13 ) 



OJ- 



for X^ f = 2{l-a), X 2 f - 



X^ q = -(1-a), i. e., for 



2.466 x 10 4 (1-a) 
TYW.& - 157.Sa)Y,{l-; 



f- V 



( 87 ) 



since 



W= 2.016a •- 159.8 (1-a) = 159.8 - 157.8a. 

The values of 7: calculated from Eqs. (86) and (87) are plotted in 

r 

Fin. 5. 



i 



* 
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VII. COMPARISON OF CALCULATED AND OBSERVED BURNING 
VELOCITIES 

In order to calculate the laminar burning velocity as a func- 
tion of mixture ratio, it is convenient to summarize first the rele- 
vant theoretical relations. 

The laminar burning velocity u^s given in terms of the param- 
eter JL through the basic relation expressed in Eq. (21), viz. , 

RT. /a7“® — 6, 172" 

T 



- .Ill 

p 

K o 



■„ nr { b i/2' p m ■ 

- — — b 2 (tt> <w> 

/a(c 17) 21j 5 K -f 



(88) 



Introduction of appropriate numerical values for E and R 

for p = 1 atmos. leads to the relation 



, „ , n 6 

u = S. 31x10 / — — — 

° / c % 



Q cm/sec 

W ° 



(89) 



if \V is expressed in g/mole, A ^ in cal/cm sec K, T q in °K, c^ in 
cal/g °K, and T, in °K. 

~ 9 2 

For the eigenvalue Ae ‘ we have obtained in Sections III, 

IV, and V the relations listed below. 

^ ! 

Bromine -Rich Mixtures (a < ■£, ^ = 2a): 



'®7 

-/Le ^ d 



(l.g }Q« f5. 1 5a-2. 10 , 2. 10(1 -a) 1 

° 2 l /TTSa m I 

— — ITJ— TTST-; T“ 

2a [ ^ 1 



(43) 



(1-1. 762a) 
(l-2a) 3/2 



< 

A e -°2 - I' 2 


0.239 VO. 761 ^ * -TZHS) — 

(1-O o ) 


d 


(ft 1.3 t 1*3*5 n 

f 'f 23 ? (»y*‘"j| 


3/2 

> 


s«', 0-« ) J 



.... ] 



( 59 ) 
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~ y 6^(1-6 Q ) [ 0. 761 + 0.7179^(l-8 o )J . 



(59a) 



Hydrogen-Rich Mixtures for 2a = 1 -a/0. 381, X. , = 2(1 -a): 

v * 



-e; 2. 32 (i-e o ) 0 J, 

vie 2 ~ — — l.V r'3. s 



f\ 1.3 , 1*30 i 



(62) 



or, in somewhat better approximation. 



-0 



Ae 



2 



1.16 

. f. T N — FT— mr- yr 

1 T ttt’TTI W"U 



(2 exr 



172 



(72) 



Hydrogen-Rich Mixtures ^ = 2-2a): 



Ae 



e'd-e ) 

-°2 °- 539 irr^r- 

~ PT Ho 



1-1. 762a 

rT72 coa 



f 1 "f£5p + (2e^p"' * ’] °* 381<a 'i ) 

* 2 






(76) 



+2 )^TH 



Explicit relations for u^can be obtained readily by combining 
Eq. (89) with Sqs. (43), (59). (72) and (76). 

Laminar burning velocities as a function of mixture compo- 
sition, computed by use of the preceding relations, are plotted in 
Fig. 6 and recorded ih Table VI. The :'.ata shown in Fig. 6 and Table 
VI correspond to the two extreme sets of values of thermal conduc- 
tivity discussed in Section VI. E. Also plotted in Fig. 6 are the ex- 
perimentally determined values for the laminar naming velocity 
reported by Anderson and his collaborators. 
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Reference to Fig. 6 shows that the calculated values for the 
burning velocity are appreciably larger than the observed results. 
Although laminar burning velocities determined from Bunsen burner 
experiments are notably unreliable, it appears likely that the major 
source of the discrepancy between theory and experiment is the re- 
sult of uncertainties in the chemical kinetics data and in the thermal 
conductivities used in our calculations. For example, on the basis 
of the ratios of rate constants listed in Table II, it is not unreason- 
able to reduce the largest calculated values of the burning velocities 
by as much as a factor of 1.3, thereby producing much better agree- 
ment between theory and experiment than is indicated in Fig. 6. 
Furthermore, our calculations are based on the steady- state approxi- 
mation for hydrogen and bromine atoms. We neglect, therefore, 

the influence of diffusion. It is well-known that the effect of diffu- 

(91 

sion on burning velocity is to decrease the burning velocity. ' ' 

Finally, it has been demonstrated that the steady-state approxima- 

J9t 

tion for bromine atoms cannot apply exactly at elevated temperatures 
for the bromine -rich and near -stoichiometric mixture compositions. 
Details concerning the corrections required because of failure of 
the steady-state approximation for Br- atoms have been worked out 
by von Karman and Millan and will not be given here. 

In spite of the discrepancies noted for ti e absolute values of 
calculated and observed burning velocities, it is evident that the cal- 
culated burning velocities have the correct functional dependence on 
composition. This result gives us confidence that, except for the 
questionable approximation concerning the steady-state approxima- 
tion for Br- atoms, we have a valid first approximation to the physico- 
chemical processes occurring in a hydrogen-bromine flame. 
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Table I. Logarithm of the equilibrium constants K ^ , for the 
reaction Br -. 2Br, as a function of temper- 
ature, talten from reference 13 

T(°K) -log K (p in atmos. ) 

P 

298.16 28.318 



300 


28. 110 


400 


19.7096 


500 


14.6528 


600 


11.2754 


700 


8.8570 


800 


7.0392 


900 


5.6220 


1000 


4.4858 


1100 


3. 5536 


1200 


2. 7746 


1300 


2.1136 


1400 


1.5454 


1500 


1.0512 


1600 


. 5730 



1700 



0175 



Table II. The ratio of reaction rates, k ., , calculated from 
gs. (79) and (79a) ae a function of temperature 



T (°K) 

500 

600 

1000 

1500 



k2 from Eg. (79)/k2 from Eg. (79a) 
.69 
.91 
1.42 
1.67 



1700 



1.72 



-43- 



Table HI. Adiabatic flame temperatures and equilibrium 
mole fractions ao a function of composition (p = 1 atmoa. 
T 0 = 323 °K, a = initial mole fraction of H£) 



1 -a 



T f (°K) 



X 



HBr, f 



X 



H 



X 



2,f 



Br 



2, f 



a, f 



k 3r,f 



. 1 675 

.2 1014 

.3 1324 

.4 1585 

.5 1683 

. 57 1475 

.6 1395 

.7 1182 

.8 923 

725 



20 


.80 


— 


— 


— 


40 


.60 


8.0x10"^ 


— 


— 


60 


.40 


— 


— 


— 


79 


.21 


1.6x1 O' 4 


2. 5xl0~” 


7 1 3x10 


95 


2.7xl0" 2 


2.9xl0" 3 


2.2xl0" 5 


4.9x10 


86 


3. 5x1 0" 4 


.085 


3.3xl0" 7 


7.9x10 


80 


7. 6xl0” 5 


. 162 


4. 3x1 0” 8 


6. 0x10 


60 


7. IxlC' 6 


.378 


1.2x1 0' 9 


3. 4x10 


40 


— 


.60 


— 


— 


20 


« • « • 


.80 


m _ ^ — 


mm mm ^ mm 



9 



Table IV. Thermal conductivities A f as a function of 



composition, calculated from the empirical relation given in 


Eq. (82) and based on the data of Campbell and Hirochfelder 1 


l-a 


A£ jcIO*''* (cal /cm sec °K) 


0 . 1 


43.2 


0.2 


39.6 


0.3 


33. 1 


0.4 


25.0 


0.5 


14. 1 


0.6 


8.8 


0.7 


6.0 


0.8 


3.6 


0.9 


2.3 



